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interplanetary  conditions  This  approach  has  the  potential  for  significantly  improving  the  performance  of  operational  models  used  to  calculate  trajectories  of 
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[i]  Orbit-averaged  mass  densities  p  and  exospheric  temperatures  T ^  inferred  from 
measurements  by  accelerometers  on  the  Gravity  Recovery  and  Climate  Experiment 
(GRACE)  satellites  are  used  to  investigate  global  energy  Eth  and  power  Ilth  inputs  to  the 
thermosphere  during  two  complex  magnetic  storms.  Measurements  show  p,  7^,  and 
E)h  rising  from  and  returning  to  prevailing  baselines  as  the  magnetospherie  electric  field 
s\ s  and  the  Dst  index  wax  and  wane.  Observed  responses  of  Eth  and  T ^  to  £Vs  driving 
suggest  that  the  storm  time  thermosphere  evolves  as  a  driven-but-dissipative 
thermodynamic  system,  described  by  a  first-order  differential  equation  that  is  identical 
in  form  to  that  governing  the  behavior  of  Dst.  Coupling  and  relaxation  coefficients  of  the 
Eih*  Toe*  and  Dst  equations  are  established  empirically.  Numerical  solutions  of  the 
equations  for  and  Eth  are  shown  to  agree  with  GRACE  data  during  large  magnetic 
storms.  Since  T ^  and  Dst  have  the  same  £Vs  driver,  it  is  possible  to  combine  their 
governing  equations  to  obtain  estimates  of  storm  time  thermospheric  parameters,  even 
when  lacking  information  about  interplanetary  conditions.  This  approach  has  the  potential 
for  significantly  improving  the  performance  of  operational  models  used  to  calculate 
trajectories  of  satellites  and  space  debris  and  is  also  useful  for  developing  forensic 
reconstructions  of  past  magnetic  storms.  The  essential  correctness  of  the  approach  is 
supported  by  agreement  between  thermosphene  power  inputs  calculated  from  both 
GRACE-based  estimates  of  Elh  and  the  Weimer  Poynting  tlux  model  originally  derived 
from  electric  and  magnetic  field  measurements  acquired  by  the  Dynamics  Explorer  2 
satellite. 
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1.  Introduction 

[2]  The  LIS  Space  Surveillance  Network  (SSN)  is  respon¬ 
sible  for  tracking  and  predicting  the  trajectories  of 
thousands  of  space  objects.  Precise  orbit  determinations 
require  accurate  knowledge  of  existing  distributions  of  mass 
in  the  thermosphere.  The  performance  of  operational  mod¬ 
els  of  the  thermosphere  needed  to  calculate  drag  forces  is 
satisfactory  during  quiet  times  but  degrades  as  geomagnetic 
activity  increases.  During  the  magnetic  storm  of  March 
1989  more  than  1500  objects  in  the  SSN  catalog  were 
temporarily  lost.  The  destruction  of  Feng  Yun  IC  by  an 
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antisatelhte  weapon  in  January  2007  created  more  than 
2500  ncwr  debris  fragments  in  low  Earth  orbit  (LEO).  The 
potential  for  debris  from  Feng  Yun  and  preexisting  objects 
to  inflict  severe  to  catastrophic  damage  on  spacecraft 
operating  in  LEO  has  refocused  attention  on  the  need  to 
strengthen  collision-avoidance  capabilities  [ Wright,  2007] 
This  in  turn  requires  improved  specifications  of  thermo¬ 
spheric  conditions  during  severe  space  weather  events. 

[3]  This  paper  describes  a  new  approach  to  specifying 
thermospheric  dynamics  that  was  developed  at  the  Air 
Force  Research  Laboratory  (AFRL)  to  support  operational 
storm  time  drag  calculations  [Bowman  el  al .,  2008].  It 
extends  three  AFRL  studies  of  ionosphere-thermosphere 
(IT)  responses  to  external  driving  during  magnetic  storms. 
The  first  study,  by  Huang  and  Burke  [2004],  showed  that 
during  the  main  phases  of  large  magnetic  storms  Defense 
Meteorological  Satellite  Program  (DMSP)  spacecraft  en¬ 
counter  repeated  episodes  of  intense  (>1000  nT)  magnetic 
perturbations  at  magnetic  latitudes  <60°.  These  events 
produce  only  weak  (<100  nT)  disturbances  on  the  ground 
Conservative  estimates  of  the  total  electromagnetic  energy 
deposited  in  the  ionosphere  over  a  typical  20-minutc  cpi- 
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sode  were  a  few  percent  of  the  ring  current  energy  deter¬ 
mined  from  the  Dessler-Parker-Sckopke  relation  [Carovil- 
lano  and  Siscoe ,  1973;  Stern ,  2005]  Huang  and  Burke 
[2004]  concluded  that:  (1)  models  that  rely  on  data  from 
ground  magnetometers  at  auroral  latitudes  greatly  underes¬ 
timate  IT  energy  inputs,  and  (2)  the  ring  current  acts  as  an 
energy  reservoir  for  heating  the  storm  time  thermosphere. 

[4]  The  second  AFRL  study  compared  these  conclusions 
with  the  dynamics  of  neutral  densities  inferred  from  meas¬ 
urements  of  accelerometers  on  the  Gravity  Recovery  and 
Climate  Experiment  (GRACE)  satellites  during  the  storms 
of  November  2004  [Burke  et  al 2007a].  This  purely 
observational  study  yielded  five  empirical  conclusions: 

[5]  1 .  Consistent  with  conclusion  ( 1 )  of  Huang  and  Burke 
[2004],  present  drag  models  use  the  ap  index  as  a  driver  and 
underestimate  observed  storm  time  density  increases  by 
>100%. 

[6]  2.  While  crossing  the  polar  cap  near  the  noon-mid- 
night  meridian  GRACE  experienced  unpredieted  large  pos¬ 
itive  and  negative  drag  spikes  as  the  spacecraft  encountered 
head  and  tail  winds  driven  by  antisunward  convecting  10ns. 

[7]  3.  Local  densities  p  vary  widely  but  orbital  averages  p 
evolve  systematically  during  storms. 

[8]  4.  Magnetospherie  elcetric  fields  £VS  ~  ^pc^LyRf 
predicted  from  parameters  observed  at  Li  anticipate  varia¬ 
tions  of  p  with  ~5  hr  lead  times.  The  symbols  4>Pc  and 
2LYRt  represent  the  polar  cap  potential  and  width  of  the 
magnetosphere,  respectively.  LY  is  the  distance  in  Rc  from 
the  Earth’s  center  to  the  magnetopause  along  the  dawn  or 
dusk  meridian. 

[9]  5.  During  main  and  early  recovery  phases  of  storms 
the  Dst  index  tracked  p  variations.  From  storm  to  storm, 
when  £*vs  — *  0,  p  relaxed  to  prestorm  levels  at  about  the 
same  rate.  Dst  returned  to  quiet  time  values  at  a  much 
slower  pace. 

[10]  Since  the  Dst  and  ap  indices  mostly  reflect  ground 
effects  of  magnetospherie  and  ionospheric  Hall  currents, 
respectively,  and  magnetospherie  electric  fields  energize 
ring  current  particles,  conclusions  (1),  (4),  and  (5)  seemed 
to  confirm  the  essential  correctness  of  the  ring  current 
reservoir  hypothesis.  Allowing  for  local  "'density  spikes” 
represented  by  (2),  the  systematic  variations  of  p  confirmed 
a  suggestion  of  Wilson  et  al.  [2006]  that,  on  a  global  scale, 
the  storm  time  thermosphere  evolves  as  a  large  thermody¬ 
namic  system  that  never  strays  far  from  diffusive  equilibrium. 

[11]  Wilson  et  al.  [2006]  analyzed  neutral  densities  in¬ 
ferred  via  the  High  Accuracy  Satellite  Drag  Model 
(HASDM)  [Casa l i  and  Barker ,  2002]  during  all  geomag¬ 
netic  disturbances  in  the  first  six  months  of  2001.  HASDM 
assimilates  drag  measurements  from  75  calibration  satellites 
with  perigees  ranging  in  altitude  from  200  to  800  km. 
Consistent  with  the  first  law  of  thermodynamics,  storm 
time  variations  in  p  reflect  the  thermal  and  gravitational 
energy  histories  of  the  thermosphere.  Compiled  energy 
histories  derived  from  HASDM  outputs  were  shown  to  be 
in  substantial  agreement  with  power  into  the  high-latitude 
ionosphere  predicted  by  the  Poynting  flux  model  of 
Weimer  [2005]  (W5).  The  W5  model  combines  statisti¬ 
cally  determined,  high-latitude  electrostatic  and  magnetic 
Euler  potential  responses  to  interplanetary  changes  to 


estimate  Poynting  fluxes  into  the  global  ionosphere  and 
thermosphere.  As  such  it  does  not  require  prior  knowl¬ 
edge  about  distributions  of  ionospheric  conductance. 
Beeause  some  fraction  of  the  energy  deposition  occurs 
at  altitudes  below  200  km,  W5  predictions  exceeded 
HASDM-based  calculations. 

[12]  In  the  third  AFRL  study  Burke  [2008]  explored 
relationships  implicit  within  the  model  of  Jacchia  [1977] 
(J77)  to  quantify  thermospheric  energy  budgets  during  large 
magnetic  storms.  The  choice  of  J77  was  one  of  conve¬ 
nience;  in  principle  these  techniques  can  be  adapted  to  any 
model  The  J77  model  uses  analytic  diffusion  formulas  to 
calculate  density  (p)  and  temperature  (T)  profiles  for  any 
exospheric  temperature  (T^).  Appendix  A  provides  the  J77 
equations  used  to  specify  temperature  profiles.  Exospheric 
temperatures  T^  uniquely  specify  both  density  and  temper¬ 
ature  profiles.  The  modeled  atmosphere  consists  of  a  fixed 
number  of  diatomic  (N2,  02)  and  monatomic  (O,  Ar,  He,  H) 
species  (ct).  At  altitudes  <90  km,  species  are  well  mixed 
and  maintain  approximately  the  ground  fractional  densities. 
A  minimum  temperature  of  188  K  is  assigned  at  the 
mesopause  altitude  h  =  90  km.  Above  90  km  all  species 
are  in  diffusive  equilibrium  specified  by  T(h)  profiles  that 
pass  through  inflection  points  at  125  km  and  approach 
at  high  altitudes.  For  energy  calculations  below  we  define, 
number  densities  n(r)  =  V  n^r),  mean  masses  m(r)  = 

nifjtia  (r),  mass  densities  ~p(r)  =  Tri(r)n(r)  and  the  heat 
capacity 

^  /  v  kBA 

•  |^(n[N2]  +  n[02])  +|(n[0]  +  n[Ar]  +  n[He]  +n[ll])j 

(I) 

where  r  =  RE  +  h  is  distance  from  the  Earth’s  center;  kB  and 
A  are  Boltzmann’s  constant  and  Avogadro’s  number, 
respectively.  The  approximation  sign  in  equation  (1) 
recognizes  the  absence  of  chemical  reactions  in  the  J77 
thermosphere.  In  the  model,  height  profiles  of  p  and  T  are 
uniquely  specified  by  T^.  Since  ~p  is  measured  and  the 
orbit-averaged  altitude  ( h )  of  GRACE  is  known,  it  seemed 
possible  to  use  these  two  parameters  to  determine  T^  Least 
squares  testing  revealed  relationships  between  T^>  and 
found  in  J77  tables  that  are  well  represented  in  the  form  of  a 
quadratic  polynomial 

2 

T^  =  Y^ai(h)f/(h)  (2) 

1=0 

In  the  300  <  h  <  500  km  altitude  range  sampled  by  the 
GRACE  and  Challenging  Minisatellite  Payload  satellites, 
the  coefficients  a^h)  are  described  by  5th-order  polynomials. 


=YLb'ih’  <3) 

/= 0 
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Specifically 


/ 1  \ 
h 


/ao(h)\ 

/  -28.10 

2.69 

-2.03  x  10  3 

0 

0 

0  \ 

a.0)  = 

-4.733  x  1017 

4.312  x  1015 

-1.372  x  1013 

1  60  x  1010 

0 

0 

V  32 [h)  / 

\  3.2695  x  1032 

-4.620  x  1030 

2.618  x  1028 

-7  456  x  1025 

1.071  x  1023 

-6  237  x  lO19/ 

Regression  coctTicients  obtained  in  fitting  T^  and  a,  to 
polynomials  exceed  0.999.  In  applying  equations  (2),  (3)  and 
(4),  p  and  h  are  in  grams/ee  and  kilometers,  respectively 
[  13]  The  total  thermospheric  energy  density  is  the  sum  of 
the  thermal  tjj{r)  =  Cy  ( r)n(r)T(r)!A ,  and  potential  energy 
contnbutions  (po  (r)  =  p(r)MfGlr  [Wilson  et  ai ,  2006].  The 
symbols  \fE  and  G  represent  the  Earth's  mass  and  the 
gravitational  constant,  respectively.  To  estimate  the  thermo¬ 
sphere's  total  energy  content  Eth  =  HT  +  <£<7,  it  is  necessary 
to  integrate  rfr(r)  and  (r)  over  the  volume  of  the 
thermosphere.  Our  use  of  p  has  effectively  averaged  latitu¬ 
dinal  variations  of  the  integrands.  Local  time  integrations 
are  effected  via  considerations  of  exospheric  temperature 
distributions  in  the  Jaeehia  models  where  T^  ranges  be¬ 
tween  some  minimum  and  maximum  values  with  a  ratio 
of  R  %  1.31.  Polar-orbiting  spacecraft  like  GRACE  sample 
globally  averaged  exospheric  temperatures  FiX,  «  1.155 7"^ 
Inin  regardless  of  the  local  times  of  their  orbital  planes.  Note 
that  current  drag  models  determine  the  parameter  Fv  min 
from  combinations  of  average  and  present  levels  of  solar 
EUV  radiance  proxies  and  the  prevailing  ap  index  [Bowman 
et  al .,  2006]  The  total  thermal  energy  is  approximately 


Kf+um  fljr+iooo 

J  Cv(r)n(r)T(r)rdr  (5) 

Rf.  +/»o  Rf.  +  ha 

Since  we  are  only  interested  in  changes  in  potential  energy 
it  is  useful  to  set  the  potential  energy  of  the  thermosphere  to 
zero  at  the  base  of  the  integration  range  [  Wilson  et  al ,  2006] 
and  represent  the  gravitational  energy  of  thermospheric 
neutrals  as 


Rt  f 1000 

<f>0  %  4 7T  I  [0(r)  -  0(ro)]r*dr 

Rf 

#*•+1000 

-  4xMt  G  J  f>(r) 

Rf 

With  information  available  in  J77  tables,  Burke  [2008] 
calculated  the  integrands  of  (5)  and  (6)  for  700  <  T^  < 
2000  K  at  increments  of  100  K  then  numerically  integrated 
them  upward  from  an  altitude  of  100  km  to  obtain  E!h  as  a 
function  of  T-*.  A  linear  regression  analysis  of  the  resulting 
Eth  as  a  function  of  F^  showed  that 


rdr 


(6) 


Etti (h  >  100km)  =  5.365  10n  4-  8.727  1013T^.  (7) 


with  a  regression  coefficient  R  >  0.998.  We  sec  that  if  F^ 
changes  by  100  K,  E,h  gains/loses  ^8.7  •  1 0 1  J.  Over  the 
range  700  <  T^  <  2000  K,  E,h  only  varies  by  ~I5%. 

[u]  The  following  section  briefly  describes  data  sources 
used  in  this  report.  The  third  section  presents  orbital 
averaged  densities  measured  by  GRACE  during  the  summer 
and  late  autumn  of  2004,  as  well  as  exospheric  temperatures 
and  thermospheric  energies  derived  from  them  via  J77. 
These  parameters  are  compared  with  Dst  and  5Vs«  Our 
analysis  of  GRACE  data  indicates  that  storm  time  E,h 
relaxes  exponentially  to  predisturbance  levels  when  £\s 
turns  off.  In  the  fourth  section  we  treat  Elh,  T^,  and  Dst 
as  manifestations  of  driven-dissipative  systems  with  the 
same  fVs  driver  and  compare  predictions  with  available 
observations.  The  fifth  section  contains  a  summary  of 
research  presented  in  this  paper,  a  list  of  its  new  conclu¬ 
sions,  and  offers  comments  on  the  limits  of  applicability  for 
coupling  coefficients  derived  from  GRACE  measurements. 
The  document  also  contains  two  appendices  that  describe 
the  J77  model  and  the  derivation  of  sv s- 

2.  Data  Sources 

[15]  Thermospheric  mass  densities  are  estimated  from 
measurements  of  nearly  identical  Spatial  friaxial  Acceler¬ 
ometer  for  Research  (STAR)  sensors  on  CHAMP  and  two 
GRACE  satellites  that  fly  in  tandem  in  nearly  circular,  polar 
orbits  [Tapley  et  al .,  2004].  STAR  sensors  monitor  electro¬ 
static  forces  needed  to  maintain  proof  masses  (PM)  at  the 
center  of  cages  located  within  2  mm  of  each  spacecraft's 
center  of  mass  [Bruinsma  et  al ,  2004].  The  spacecraft  and 
PM  respond  to  gravity  in  the  same  way.  Thus  changes  in  the 
electrostatic  forces  that  maintain  the  PM  at  the  center  of  its 
cage  reflect  spacecraft  responses  to  nongravitational  forces 
such  as  atmospheric  drag  and  radiation  pressure  [Bruinsma 
and  Biancale ,  2003].  Acceleration  due  to  atmospheric  drag 
is  given  by 


adng  =  CD(A*/M>c)pV2  (8) 

where  Asc  and  Msc  represent  the  cross-sectional  area  and 
mass  of  the  spacecraft,  respectively,  p  the  mass  density  of 
the  neutral  atmosphere,  and  V  is  the  spacecraft  velocity  in 
the  rest  frame  of  ambient  neutrals.  The  drag  coefficient  Cn 
depends  on  the  angle  of  flow  to  the  spacecraft  surface,  the 
ratio  of  the  temperatures  of  the  satellite  surface  and  the  local 
atmosphere,  and  the  ratio  of  the  mean  mass  of  atoms  in  the 
atmosphere  to  those  on  the  satellite  surface  [Bruinsma  and 
Biancale ,  2003] 

[16]  To  maintain  the  same  temporal  eadenee  as  Dst,  5Vs 
values  presented  in  this  paper  were  calculated  using  hourly 
averaged  values  of  solar  w  ind  and  IMF  parameters  from  the 
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Figure  1.  Geophysical  parameters  and  thermospherie 
responses  observed  by  GRACE.  From  top  to  bottom,  plots 
show  the  history  of  the  Dst  index,  £Vs>  as  well  as  orbit- 
averaged  measurements  of  mass  density  J\  exosphenc 
temperature  T ^  and  thermospherie  energy  Eth.  Data  are 
plotted  as  functions  of  universal  time  from  Julian  days  150 
to  230,  2004. 


Advanced  Composition  Explorer  (ACE)  satellite  in  a  halo 
orbit  around  the  L!  point,  ~234  Rp  upstream  of  Earth.  The 
solar  wind  density  and  velocity  are  measured  by  the  Solar 
Wind  Electron,  Proton,  and  Alpha  Monitor  (SWEPAM) 
[Me  Co  mas  et  al .,  1998].  The  interplanetary  magnetic  field 


(IMF)  vector  is  observed  by  the  magnetic  field  instrument 
(MFI)  [Smith  et  al,  1998]. 

3.  GRACE  Observations 

[17]  Figure  1  includes  plots  of  the  geophysical  parameters 
Dst  (top)  and  £Vs  (second)  as  well  as  orbit-averaged 
measurements  of  p  (third),  T ^  (fourth),  and  Eth  (bottom) 
acquired  by  GRACE  between  1  June  and  17  August  2004 
(Julian  days  (JD)  150  to  230).  During  this  interval 
GRACE’S  orbit-averaged  altitude  was  ~486.5  km;  its 
orbital  plane  precessed  from  near  the  noon-midnight  to  the 
dawn-dusk  meridian.  Values  of  T ^  and  E^  were  obtained  by 
applying  equations  (2)  and  (7)  to  the  p(t)  stream.  The  data 
plots  show  that  T x  ranged  from  750  to  1200  K  and  Erh 
between  6.05  and  6.4  •  1017  J  over  this  period. 

[18]  Last  three  data  plots  in  Figure  1  show  that  p(ty  7\ 
and  Efh  vary  on  two  time  scales:  (1)  low-frequency  varia¬ 
tions  with  the  27-day  solar  rotation,  and  (2)  outerops  that 
emerge  in  response  to  changes  in  geophysical  conditions 
related  to  intensifications  of  Dst  and  eVs-  The  largest  change 
in  the  Eth  trace  of  ~3*  10 16  J  occurred  on  JD  209  in  response 
to  interplanetary  forcing  when  Dst  attained  a  minimum 
value  of  about  -200  nT.  The  energy  of  the  ring  current 
Erc  estimated  via  the  Dcsslcr-Parker-Sckopke  (DPS)  rela¬ 
tion  [Stern y  2005]  is 

Erc(J)  *  W  2.58  •  1013  Dst(nT)  (9) 

2  no 

where  B0  ~  3.1  •  104  nT  is  the  magnetic  field  at  the  equator, 
and  Em  ~  8  •  10 1  J  is  the  energy  of  the  magnetic  field 
above  the  Earth’s  surface.  Equation  (9)  indicates  that  at  its 
maximum  value  ERr  «  5.16  •  101*  J  and  was  substantially 
smaller  than  the  increase  in  E,h.  This  empirieal  result 
indicates  that  the  ring  current  reservoir  hypothesis  is 
untenable.  A  large  fraction  of  the  storm  time  energy  budget 
comes  from  the  solar  wind  to  the  thermosphere,  unmediated 
by  the  ring  current. 

[ 1 9]  Figure  2  compares  the  same  geophysieal  and 
GRACE  parameters  measured  during  JD  305-330  (2-27 
November  2004)  using  the  format  of  Figure  l .  The  orbital 
plane  of  GRACE  was  near  the  noon-midnight  meridian.  The 
Dst  plot  shows  that  two  major  geomagnetie  disturbances 
oeeurred  between  JD  312  and  315  (7-10  November). 
Similar  disturbances  appear  on  all  of  the  other  plots.  Dst 
minima  were  -373  nT  (06:00  UT,  8  November)  and 
-289  nT  (10:00  UT,  10  November).  The  corresponding 
values  of  ER£  are  9.62  and  7.46  •  10IS  J;  predisturbanee 
baselines  for^^  and  Eth  were  near  900  K  and  6.15  •  101  J, 
respectively.  rose  to  1390  K  at  07:00  UT  on  JD  3 13  and 
1350  K  at  10:00  UT  on  JD  315.  Eth  maxima  at  these  times 
were  6.58  and  6.52  •  101  J  Again,  corresponding  AEth 
values  of  4.3  and  3.7  •  1016  J  exceed  estimates  of  ERr  9.6 
and  7.5  •  10 15  J,  respectively. 

[20]  Burke  et  al.  [2007a]  found  that  storm  time  thermo¬ 
spherie  density  increases  always  relaxed  at  the  same  rate 
when  £Vs  turned  off.  It  turns  out  that  the  storm  time 
thermosphere  loses  its  storm  time  energy  increment  expo¬ 
nentially  with  time.  This  is  demonstrated  in  Figure  3  in 
which  £Vs  and  the  natural  logarithm  of  E,h  are  plotted  as 
functions  of  UT  during  the  disturbed  interval  JD  204  -21 1 
(22-29  July  2004).  Vertical  lines  mark  times  when  £Vs 
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Figure  2.  Geophysical  parameters  and  thermospheric 
responses  observed  by  GRACE  plotted  as  functions  of 
time  for  JD  306  to  330,  2004  using  the  same  format  as 
Figure  1 . 


rapidly  decreased  to  low  values.  The  two  slanted  lines  that 
overlay  the  trace  of  In  Eth  during  the  subsequent  periods  of 
low  have  identical  slopes.  Numerically,  Eth  has  an 
enfold  relaxation  time  r  %  6.5  hours. 


4.  Modeling  a  Driven-Dissipative  Thermosphere 

[21]  The  behavior  of  Elh,  responding  to  variations  in  eVs 
and  decaying  exponentially  to  predisturbance  levels  when  it 


turns  off,  is  reminiscent  of  commonly  observed  phenomena 
called  driven-dissipative  systems.  The  magnetic  energy 
generated  by  a  current  flowing  in  a  resistor-inductor  (R-L) 
circuit  provides  a  familiar  example  [Young,  1992].  Burton  et 
al.  [1975]  proposed  that  the  ring  current  evolves  in  a  similar 
way  and  derived  a  linear  differential  that  describes  changes 
in  Dst*,  the  pressure  corrected  Dst. 


dDst*  Dst* 

- T—  =  I - 

at  r  RC 


(10) 


Here  5/,  aD,  and  trc  represent  the  interplanetary  electric 
field,  the  coupling  coefficient  and  the  relaxation  time 
constant  of  the  ring  current,  respectively. 

[22]  The  following  subsections  consider  the  driven-dissi¬ 
pative  properties  of  Elh,  T^  and  Dst.  The  analogy  of  rolling 
hills  and  rocky  outcrops  drawn  from  GRACE  data  in  Figure 
1  indicates  that  we  may  regard  the  scalar  fields  Eth  and  T x 
as  having  two  independent  additive  sources,  the  solar  EUV 
radiance  and  the  solar  wind.  The  dnven-dissipative  equa¬ 
tions  presented  below  describe  the  solar  wind  contributions 
to  these  scalars. 

4.1.  Energy  and  Power 

[23]  During  magnetic  disturbances  E,h  variations  correlate 
with  £vs-  Rather  than  pursue  statistical  relationships  be¬ 
tween  these  two  quantities,  we  propose  to  exploit  implica¬ 
tions  of  information  contained  in  the  decay  characteristics  of 
Eth  illustrated  in  Figure  3.  Specifically  we  postulate  that  on 
a  global  scale  the  storm  time  thermosphere  acts  like  a 
driven-dissipative  system.  If  this  conjecture  is  correct,  the 
solar  wind  contribution  to  Eth  should  evolve  in  ways  similar 
to  Dst  in  equation  (10).  The  governing  equation  takes  the 
form 


dEthsw 

dt 


—  (*e£vs  “ 


Eihsw 
T E 


(11) 


Figure  3.  Plots  of  eVs  and  the  natural  logarithm  of  Eth  sw 
for  the  disturbance  on  JD  206  to  211,  2004.  Vertical  lines 
mark  times  of  electric  field  decrease.  The  slanted  lines  have 
the  same  slopes  indicating  that  Eth  sw  decays  exponentially 
when  £Vs  turns  off.  The  estimated  e-fold  relaxation  time  r  is 
M).5  hours. 
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thermospherie  energy  changes  with  independent  W5  pre¬ 
dictions  of  electromagnetic  power  into  the  thermosphere. 
W5  uses  solar  wind/IMF  measurements  to  predict  distribu¬ 
tions  of  Poynting  flux  into  the  northern  and  southern  iono¬ 
spheres.  Integration  over  the  affected  areas  provides  the 
global  rates  of  energy  input  PW5.  On  the  other  hand,  orbit- 
averaged  Eth  measurements  from  GRACE  are  snapshots  that 
capture  the  thermosphere’s  total  energy  content  averaged 
over  orbital  periods.  GRACE  data  reflect  balances  that  the 
thermosphere  strikes  between  power  received  from  the  solar 
wind  and  lost  via  radiative  and  other  processes  [Mlynczak  et 
al ,  2005]. 

[26]  Equation  (11)  provides  a  bridge  between  the  two 
approaches.  The  term  apsys  represents  the  rate  at  which 
energy  is  provided  to  drive  the  storm  time  thermosphere. 
Thus  it  can  be  compared  with  PW5  directly.  Since  aF  ~  5.5  • 
I015  [(J/hr)/(mV/m)]  and  P\v5  is  in  Watts  wc  must  divide  aF 
by  3600  s/hr  to  obtain  qe  %  1.528  •  1012  W/(mV/m)  or 
1.528  TW/(mV/m).  Figure  5  shows  plots  of  PW5  (red)  at  a 
5-minute  cadence  and  hourly  averaged  II th  «  aF  £vs  in  TW 
for  the  storms  of  July  and  November  2004.  The  qe  £vs  11306 
was  time  shifted  by  1  hour  to  allow  for  transport  from  ACE 
to  the  magnetosphere.  Delay  times  used  in  W5  vary 
according  to  the  speed  of  the  solar  wind  and  the  tilt  angle 
of  surfaces  of  constant  interplanetary  electric  field  phase 
\Weimer  et  al .,  2002].  In  general  the  agreement  between 
P W5  and  the  power  required  to  explain  GRACE-based 
calculations  of  Elh  appears  remarkably  good.  Thus,  on  a 


The  coupling  constant  aF  and  relaxation  time  rF  =  6.5  hours 
are  empirically  determined.  Equation  (7)  can  be  solved 
numerically 

EthS\v(tnf  1 )  =  0£-5Ks(tn)At  +  Ethsw(tn)  ^  1  “  (12) 

where  At  represents  the  time  step  between  samples.  For 
convenience  we  let  At  =  1  hour,  and  express  £ys  in  mV/m. 
In  this  ease 


Ethsw(tn+l)  —  £U*£Vs(tn)  +  £thsw(tn)(  1  “  —j 

=  ^£^vs(tn)  +  .846EthSw(tn)  (13) 

Through  trial  and  error  comparisons  with  GRACE  measure¬ 
ments  acquired  during  JD  150-230,  2004  we  found  that  aF 
^  5.5  •  101S  [(J/hr)/(mV/m)].  Numerical  solutions  of 
equation  (13)  can  be  compared  with  GRACE  measurements 
to  confirm  or  refute  the  assumption. 

[24]  Figure  4  contains  plots  of  £vs  (black),  solutions  of 
equation  (13)  (blue)  and  GRACE  measurements  of  Eth  sw 
(red).  The  plots  represent  the  disturbed  periods  in  (top)  July 
and  (bottom)  November  2004.  In  both  instances  we  approx¬ 
imated  Eth  sw  as  Eth  —  6. 1  •  10l  J  since  the  E^  sw  calcula¬ 
tions  in  Figure  5  are  presented  in  units  of  1016  J.  Although 
aE  was  determined  using  data  from  the  summer  of  2004 
only,  predictions  of  equation  (13)  also  appear  to  be  in  good 
agreement  with  GRACE  measurements  during  November 
2004. 

[25]  The  essential  correctness  of  equation  (13)  can  also  be 
verified  or  falsified  by  comparing  GRACE  observations  of 
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Figure  5.  Comparison  of  storm  time  power  into  the  global 
thermosphere  predicted  by  the  W5  model  (red)  and  as Vs 
(black)  plotted  as  functions  of  UT  during  magnetically 
disturbed  periods  in  (top)  July  and  (bottom)  November 
2004. 
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Figure  4.  Comparison  of  svs  (black),  modeled  Eth  sw 
(blue),  and  Eth  sw  values  inferred  from  GRACE  measure¬ 
ments  of  ~p  (red  dots)  plotted  as  functions  of  universal  time 
during  magnetically  disturbed  periods  of  (top)  July  and 
(bottom)  November  2004.  E^  sw  is  plotted  in  units  of 
1016  J. 
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Figure  6.  Model  predictions  of  Tl3c,  $\V  (black)  and  fvs 
(red)  inferred  from  hourly  averaged  parameters  from  ACE, 
plotted  as  functions  of  universal  time  for  JD  150  to  230, 
2004. 


global  scale,  power  inputs  predicted  by  W5  and  needed  to 
model  GRACE  measurements  are  viewed  as  mutually 
consistent.  Consequently,  PW5  as  well  as  aE  £Vs  can  be 
used  as  a  driver  in  equation  (11)  to  estimate  the  develop¬ 
ment  of  E,h  sw- 

[27]  Attention  is  directed  to  two  particular  aspects  of  the 
results  shown  in  Figure  5. 

[2k]  1.  GRACE  measurements  and  W5  predictions  indi¬ 
cate  sustained  electromagnetic  power  input  >1  TW  during 
both  the  July  and  November  storms.  This  is  significantly 
greater  than  the  storm  time  inputs  due  to  energetic  particle 
precipitation  and  from  solar  UV  [Knipp  et  a/.,  2005]. 

[29]  2.  During  the  third  disturbance  of  the  July  period  and 
the  first  one  of  the  November  storms,  PWs  >  c*e  £vs  by 
significant  amounts  This  discrepancy  probably  reflects  the 
different  estimates  of  4>PC-  In  both  instances  the  polar  cap 
potential  used  in  W5  was  larger  than  that  predicted  by  the 
Siscoc-Hill  model  (Appendix  A).  For  a  given  system  of 
field-aligned  currents  the  one  driven  by  the  higher  <I>Pt 
requires  a  higher  power  input  to  overcome  eollisional  drag 
exerted  by  neutrals  on  E  x  B  drifting  ions. 


4.2.  Exospheric  Temperatures 

[30]  While  Eth  may  be  of  interest  for  estimating  energy 
transfer  from  the  solar  wind  and  magnetosphere  to  the 
thermosphere  via  the  high-latitude  ionosphere,  T x  is  gen¬ 
erally  a  more  useful  parameter  for  operational  modelers  of 
atmospheric  drag  [Marcos  et  ai ,  1998].  The  linear  relation¬ 
ship  between  E,h  and  T x  captured  in  equation  (7)  indicates 
that  a  differential  equation  similar  to  (1 1)  also  describes  the 
evolution  of  We  assume  that  T x  has_independent 

UV  and  solar  wind  sources,  7'oc  =  T no  uv  +  T.x  sw-  Since 
the  solar  UV  provides  thermospheric  baselines 


7%cuv  - 


Eth uv(J)  -  5.365  101 
8.727-  1013 


(14) 


and 


7=  _  Ethsw(J) 

xSW  8.727  I0n 

substitution  of  (15)  into  (13)  gives 

Tv sw(tn+ 1 )  =;  or-vs(tn)  +  Tx  sw(tn)  ^1  —  — 
~  63fvs(|«)  +  .846Tx.sw(tn) 


(15) 


(16) 


The  coupling  coefficient  aT=  aEIS.121  ■  1013  =  5.5  •  1 0 1 5/ 
8.727  •  10n  «  63. 

[31]  Figure  6  shows  plots  of  £vs  (red)  and  T x  sw  (black) 
for  JD  150-230,  2004.  Attention  is  directed  to  three  aspects 
of  these  data: 

[32]  1.  Whenever  solar  wind  data  from  ACE  were  un¬ 
available  we  set  6\s  =  0-  In  these  intervals  modeled  S\y 
decreased  exponentially,  but  quickly  recovered  when  the 
£Vs  data  stream  resumed. 

[33]  2.  During  periods  of  relative  magnetic  quiet  T x  sw 
ranged  between  70  K  and  100  K. 

[34]  3.  Spikes  in  eVs  always  produced  analogous 
responses  in  sw- 

[35]  Figure  7  compares  variations  of  modeled  T ^  SVv 
(blue)  with  values  inferred  via  equation  (2)  from  measure¬ 
ments  of  p  by  accelerometers  on  GRACE  during  the  two 
magnetic  storm  periods.  For  the  July  and  November  storms 
we  estimated  sw  by  subtracting  constant  baseline 
estimates  of  900  K  from  the  full  measured  values  of  Tx 
shown  in  Figures  1  and  2.  Again,  agreement  between  model 
predictions  and  observations  appears  quite  good.  Closer 
agreement  between  data  and  the  model  can  be  achieved 
using  time-varying  representations  of  Tx  uv  baselines  that 
allow  for  solar  rotation  induced  trends  apparent  in  GRACE 
data  [Bowman  et  a/.,  2008]. 


Figure  7.  Modeled  T^sw  (blue)  and  values  inferred  from 
GRACE  measurements  of  p  (red  dots),  plotted  as  functions 
of  UT  during  the  magnetically  disturbed  periods  of  (top) 
July  and  (bottom)  November  2004.  T^Sw  was  approxi¬ 
mated  by  subtracting  850  K  from  GRACE-based  estimates 

of  Toe 
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Figure  8.  Comparison  of  (black),  modeled  (blue),  and 
measured  (red)  Dst  plotted  as  functions  of  UT  during  the 
magnetically  disturbed  periods  of  (top)  July  and  (bottom) 
November  2004.  For  ease  in  making  compansons  with  sv s> 
plots  show  the  negative  of  Dst. 


time  magnetosphere  (cf.  Appendix  B).  To  test  the  utility  of 
making  this  substitution  Figure  8  compares  £vs  (black)  with 
solutions  for  equation  (17)  (blue)  and  provisional  Dst  (red) 
as  functions  of  time  during  the  July  and  November  2004 
storms.  While  lacking  the  precision  and  sophistication  of 
artificial  intelligence  solutions  presented  by  Temerin  and  Li 
[2002,  2006]  equation  (17)  gives  an  approximate  represen¬ 
tation  of  Dsfs  storm  time  development, 

[37]  Over  the  past  decade  the  solar  wind/IMF  data  needed 
to  calculate  £Vs  have  been  available  almost  continuously. 
However,  critical  input  data  were  often  unavailable  in  the 
past  and  may  become  so  in  the  not-too-distant  future. 
During  many  storms  energetic  solar  protons  render  it 
impossible  to  determine  solar  wind  densities  in  real  time. 
Nevertheless  operational  responsibility  for  calculating  the 
precise  trajectories  of  space  objects  continues.  Knowledge 
of  is  critical  for  estimating  thermospheric  densities  and 
the  corresponding  atmospheric  drag  (equation  (8)).  The 
remainder  of  this  subsection  explores  the  possibility  of 
using  Dst  to  provide  real-time  information  regarding 
changes  in  during  magnetic  disturbances. 

[38]  It  is  possible  to  exploit  the  formal  similarities  be¬ 
tween  differential  equations  for  T^  and  Dst  to  calculate  the 
enhanced  storm  time  drag  on  space  objects  even  in  the 
absence  of  data  from  Lp  This  is  accomplished  by  eliminat¬ 
ing  the  common  term  ev s  from  the  two  equations 


«T  . 


V  (tn+ 1 ) 


Op  _ 


Dst(tn+l) 


-R) 
-(R) 


TxSw(tn) 
Dst(tn) 


=  -vs(tn) 


Rearranging  terms  gives 


(18) 


4.3.  Thermosphere-Dst  Relationships 

[36]  Burke  et  al.  [2007a]  demonstrated  that  during  mag¬ 
netic  storms  the  Dst  index  evolves  in  ways  that  closely 
approximate  the  forms  of  p(t)  traces  throughout  the  main 
and  early  recovery  phases  of  storms.  In  late  recovery  phases 
Dst  relaxes  over  several  days  while  p(t )  returns  to  its 
prestorm  baseline  in  a  few  hours.  Burke  et  al.  [2007b] 
found  that  during  magnetic  storms  with  simple  main  phases, 
Dst{t)  =  a  +  b  •  l vs  (Otwith  very  high  correlation  coef¬ 
ficients,  where  IVs  (0  =  /  £vs  (t  )df.  This  approach  ignores 

to 

the  main  phase  dissipation  term  in  equation  (10).  Linear 
regression  slopes  varied  over  a  relatively  wide  range.  The 
steepest  slopes  of  about  -40  nT/(mV/m)  were  obtained 
during  storms  with  sharp  transitions  between  main  and 
recovery  phases.  We  approximate  (10)  as 


dDst 

— 7—  =  C*£)£vs 
at 


Dst 

7*RC 


(17) 


using  aD  =  -40  [(nT/hr)/(mV/m)]  and  rRr  =  7.7  hour 
[http://sprg.ssl.berkeley.edu/dst_  index/welcome. html].  This 
approach  has  two  obvious  advantages:  (1)  Eth,  and  Dst 
can  all  be  regarded  as  responding  to  the  same  force  field, 
and  (2)  sv s  approximates  the  electric  field  that  energizes 
ring  current  particles  and  incorporates  the  essentially 
nonlinear  coupling  of  the  solar  wind/IMF  with  the  storm 


L/On  +  l) 


Tx>sw(tn) 


Dst(tn  +  I) 


(19) 


To  solve  (10)  numerically  wc  set  sw  (/0)  =  0;  the  ratio 
nj k\o  ~  1 .575.  Since  we  are  only  interested  in  the  storm 
time  Dst,  following  Burton  et  al.  [1975]  wre  set  the  '  quiet 
time”  Dst  as  —20  nT  and  impose  the  numerical  constraint 


Dst:  ifDst<-20nT 
Ust  “  —  20nT  if  Dst  >  —  20nT 


(20) 


Figure  9  plots  calculated  values  of  T c*  sw  (0  (blue) 
obtained  using  the  Dst  time  series  as  the  driver  of  equation 
(19).  The  red  dots  again  represent  values  of  7\  sw  =  T “ 
900  K  inferred  from  GRACE  measurements  of  ~p.  By  and 
large  agreement  between  the  model  and  data  appears  quite 
good.  Note  that  on  JD  210  and  316  the  modeled  T ^  sw  (/) 
decays  at  a  much  slower  rate  than  is  actually  observed.  This 
reflects  the  slowr  decay  rate  of  the  symmetric  ring  current  in 
the  late  main  phase  when  electric  coupling  between  the 
magnetosphere  and  ionosphere-thermosphere  weakens. 
Bowman  et  al.  [2008]  describe  practical  ways  to  overcome 
this  deficiency. 

[39]  It  is  also  possible  to  apply  this  technique  to  make 
forensic  reconstructions  of  thermospheric  and  interplanetary 
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Julian  Da\:  2004 


Figure  9.  Predicted  T^  Sw  (blue)  with  Dst  as  driver  and 
values  inferred  from  density  measurements  by  GRACE 
(red)  plotted  as  functions  of  UT  during  the  magnetically 
disturbed  periods  of  (top)  July  and  (bottom)  November 
2004.  As  in  Figure  7,  measured  T^  Sw  was  approximated 
by  subtracting  850  K  from  GR ACE-based  estimates  of  T^. 


conditions  during  previous  large  magnetic  disturbances 
when  solar  wind  measurements  were  unavailable.  Because 
of  its  widespread  space  weather  effects  the  March  1989 
storm  offers  an  interesting  example.  The  trace  in  the  bottom 
plot  of  Figure  10  shows  that  early  on  14  March  Dst 
approached  a  minimum  of  -600  nT.  Figure  10  (middle) 
shows  T oc  sw  (0  obtained  by  £  numerical  solution  of 
equation  (19).  We  estimate  that  T ^  rose  ~800  K  above 
prestorm  background  levels,  much  more  than  the  300  to 
400  K  increases  of  July  and  November  2004.  Figure  10 
(top)  plots  an  estimate  of  £Vs  obtained  through  a  numerical 
solution  of  equation  (18).  Near  the  maximum  epoch  of  the 
main  phase  the  calculated  £Vs  exceeded  4  mV/m.  Although 
this  value  is  large  in  comparison  with  the  July  and  November 
2004  storms,  it  is  not  outlandish.  During  the  storm  of 
March  1991  the  Combined  Release  and  Radiation  Effects 
Satellite  (CRRES)  measured  sustained  dawn-to-dusk  electric 
fields  of  mV/m  in  the  inner  magnetosphere  [H ygant  et 
ai ,  1998]. 

5.  Discussion  and  Conclusions 

[40]  This  section  summarizes  the  research  reported  above, 
lists  its  three  new  scientific  conclusions,  and  comments  on 
its  range  of  applicability. 

[41]  During  the  main  phases  of  large  magnetic  storms 
DMSP  spacecraft  regularly  cross  sheets  of  intense  (>1  A/m) 
field-aligned  currents  at  auroral  latitudes  that  produce 
very  weak  ground  signatures.  Energy  deposited  over  the 
^20-minute  lifetimes  of  these  structures  accumulates  to  a 
significant  fraction  of  the  ring  current  energy.  Huang  and 
Burke  [2004]  suggested  that:  (1)  models  using  ground 
magnetometer  measurements  underestimate  energy  inputs 
to  the  storm  time  thermosphere,  and  (2)  the  ring  current  acts 


as  an  energy  reservoir  from  which  the  thermosphere  draws 
energy.  Direct  comparisons  of  GRACE  measurements  with 
predictions  of  presently  used  models  confirmed  the  first 
conclusion  [Burke  et  ai ,  2007a].  Observed  correlations 
between  variations  of  ’p  and  Dst/£Vs  seemed  to  confirm  the 
second.  In  this  paper  we  applied  a  new  technique  to  deter¬ 
mine  exospheric  temperatures  and  thermospheric  energy 
contents  from  GRACE  measurements  of  /}.  Data  presented 
in  Figures  1  and  2  show  that  during  large  magnetic  storms  the 
energy  content  of  the  global  thermosphere  rises  by  several 
times  that  of  the  ring  current  estimated  via  the  Dessler- 
Parker-Sckopke  relation.  The  observed  responses  of  Eth  to 
Sys  led  us  to  explore  the  possibility  that  the  storm  time 
thermosphere  and  ring  current  act  as  driven-dissipative 
systems  with  the  same  driver  but  different  relaxation  time 
constants.  Numerical  solutions  of  the  differential  equation 
governing  the  behavior  of  such  systems  were  compared  with 
the  observed  variability  of  Elh,  T^,  and  Dst.  Agreement 
between  data  and  model  predictions  confirms  the  essential 
correctness  of  the  driven-dissipative  perception.  Given  their 
widely  different  data  sources,  we  regard  agreement  between 
estimates  of  electromagnetic  power  into  the  thermosphere 
predicted  by  ap£Vs  and  Pw5  as  strong,  independent  confir 
mation  of  our  dnven-dissipative  perception  of  the  storm  time 
thermosphere. 


Figure  10.  Model  estimates  of  (top)  £Vs  and  (middle)  Tv 
sw  driven  by  (bottom)  Dst  during  the  magnetic  storm  of 
March  1989. 
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[42]  The  three  major  new  conclusions  of  this  paper  are. 

[43]  1.  Storm  time  energy  increments  to  the  global 
thermosphere  are  well  represented  as  those  of  a  driven- 
dissipative  thermodynamic  system. 

[44]  2.  Rates  of  electromagnetic  energy  inputs  to  the 
thermosphere  consistent  with  p  variations  measured  by 
GRACE  match  the  independent  predictions  of  the  W5  model. 

[45]  3.  Since  T^  sw  and  Dst  respond  to  the  same  5Vs 
driver,  quantitative  estimates  of  changes  in  the  storm  time 
thermosphere  can  be  extracted  even  when  interplanetary 
data  are  unavailable. 

[4 f>]  We  close  with  a  caveat  emptor.  During  the  research 
period  leading  up  to  the  writing  of  this  paper  our  access  to 
reduced  CHAMP  and  GRACE  measurements  was  limited. 
Since  the  CHAMP  data  segment  had  significant  gaps  during 
magnetic  storms,  our  analysis  concentrated  on  measure¬ 
ments  from  GRACE  during  the  summer  and  late  fall  of 
2004.  Thus  this  report  is  more  an  extended  case  study  than  a 
statistical  analysis.  Our  estimate  of  aE  «  5.5  101<r  [(J/hr)/ 

(mV/m)]  relied  solely  on  GRACE  measurements  acquired 
during  the  July  2004  storm  period.  With  no  adjustments  this 
value  also  replicated  GRACE-based  estimates  of  Eth  sw 
during  the  November  storms.  Because  aE  £Vs  and  the 
independently  determined  Pws  are  in  substantial  agreement 
during  the  two  storm  periods  wc  regard  this  estimate  as  well 
grounded. 

[47]  That  said,  we  also  regard  the  stated  value  of  qe  as 
approximate  rather  than  definitive.  The  main  reason  for  this 
assessment  derives  from  the  recent  analysis  of  Bowman  et 
al.  [2008]  who  demonstrated  that  the  independent  calibra¬ 
tions  of  accelerometers  on  CHAMP  and  GRACE  are 
irreconcilable.  The  calibration  differences  are  not  large. 
However,  it  will  take  time  to  reach  consensus  on  the  best 
in  situ  calibrations  of  the  CHAMP  and  GRACE  acceler¬ 
ometers. 


Appendix  A:  Model  of  the  Thermosphere  of 
Jacchia  (1977| 

[4«]  J77  is  a  static  model  of  thermospheric  densities  based 
on  analytically  defined  temperature  profiles.  It  approximates 
the  upper  atmosphere  as  being  composed  of  N2,  02,  O,  He, 
and  H  with  all  species  sharing  the  same  temperature 
profiles.  These  gases  are  considered  well  mixed  below 
100  km  and  in  diffusive  equilibrium  above  it.  The  meso- 
pause  is  set  at  the  height  z*,  =  90  km  where  J77  imposes  the 
boundary  conditions:  pw  =  3.43  x  10~9  g/cm\ Tq0-  188  K, 
and  dT^dz  =  0.  At  zx  =  125  km  all  temperature  profiles 
have  inflection  points  (d2Tx/dz2  =  0),  where  dTx/dz  is 
continuous.  Above  125  km  T  increases  and  asymptotically 
approaches  an  exospheric  temperature  T*,.  Both  the  tem¬ 
perature  Tx  and  dTx/dz  at  the  inflection  point  arc  functions 
of  T^. 

[49]  In  the  J77  model  the  relationship  between  Tx  and 
is  given  by 

Tx  =  188+  WO.SSinh  *[.0045(7^  -  188)]  (Al) 
For  90  <  z  <  125  km 


T  =  Tx  4-  A  tan " 1 


For  z  >  125  km 


T=  Tt  +  Alan-'  (z  -  r.)  [l  +  5.5  10  5(=  -  rt)2]  }  (A3) 

where  Gx  =  1 .9(^3)  and  A  =  (2/ttXT^  -  Tx). 

[50]  Above  100  km,  densities  of  individual  species  n(i) 
are  described  by  the  diffusion  equation 


dn(i) 

n{i) 


M,g  /-  dTn  j.  ^ 

-^-yd+o.) 


(A4) 


where  Mj  is  the  molecular  weight  of  the  ith  species,  g  is  the 
local  acceleration  due  to  gravity,  R*  is  the  universal  gas 
constant  8.31432  joules/(mole  °K),  and  the  thermal 
diffusion  coefficient  «,  =  -0.38  for  helium  and  or;  =  0  for 
other  constituents. 


Appendix  B:  Yrolland-Stern  Electric  Fields 

[si]  As  formulated  by  Ejiri  [1978]  the  Volland-Stem  (V- 
S)  model  offers  a  simple,  albeit  limited,  method  for  esti¬ 
mating  the  trajectories  of  ring  current  ions  in  the  inner 
magnetosphere.  One  difficulty  in  applying  the  model 
concerned  the  connection  between  the  interplanetary  elec¬ 
tric  field  (IEF)  and  its  magnetospheric  manifestations.  Burke 
[2007]  reformulated  V-S  by  combining  it  with  the  Siscoe- 
Hill  (S-H)  model  of  the  polar  cap  potential  (<J>PC)  to  show 
that  in  the  absence  of  shielding  the  electric  field  in  the  inner 
magnetosphere  is  given  by  5Vs  %  <Fpc/2ReLy,  where  the 
denominator  represents  the  width  of  the  magnetosphere 
along  the  dawn-dusk  line.  This  appendix  outlines  our 
method  for  calculating  values  of  5Vs  used  in  the  main  text. 

[52]  Empirical  studies  suggest  that  on  the  dayside  the 
equatorial  magnetopause  is  nearly  self-similar  in  shape  with 
Lv  %  1  5  Lx  [Roelof  and  Si  beck,  1993].  Force  balance  at  the 
subsolar  magnetopause  requires  that 


Lx 


9.6 


llopsH  tyPSn(nPa) 
14.4/v/  Psw(nPa) 


and 


(Bl) 


More  precise  estimates  of  magnetopause  dimensions  that 
include  erosion  effects  can  be  introduced  where  more 
rigorous  calculations  are  required  [Yang  et  al 2003]. 

[53]  Between  1980  and  2000  many  investigators  sought 
linear  relationships  between  4>PC  and  interplanetary  param¬ 
eters.  Burke  et  al.  [1999]  surveyed  published  analyses  of 
4>pc  dependence  on  solar  wind  parameters  and  showed  that 
in  the  linear  regime,  consistently  high  correlation  coeffi¬ 
cients  were  obtained  with  the  relationship 


<MkV)  =  <I>0  +  LcVswBxSin2  ^  (B2) 

The  subscript  E  refers  to  IEF  contributions;  4>0  is  a  residual 
potential  with  typical  values  between  20  and  30  kV,  Bj  = 
y/(By  +  B7),  and  0  =  Cos  '(B7/BT)  is  the  IMF  clock  angle 
in  the  Y-Z  plane.  The  term  VSn  Br  Sin2  012  roughly 
corresponds  to  £\  in  equation  (10)  [Burton  et  al.,  1975].  If 
we  express  Vsw  BT  in  kV/RE  (where  1  mV/m  ^  6.4  kV/RE) 
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then  the  regression  slope  Lo  represents  a  3  to  4  RE  wide 
“gate”  in  the  solar  wind  through  whieh  geoeffeetive 
streamlines  (equipotentials)  must  pass  to  reaeh  the  dayside 
magnetopause. 

[54]  Siscoe  et  ai.  [2002]  suggested  that  during  large 
magnetic  storms  Region  1  currents  generate  perturbation 
magnetic  fields  that  alter  the  shape  of  the  dayside  magne¬ 
topause,  thereby  limiting  aeeess  of  solar  wind  streamlines 
and  forcing  the  merging  rate  to  saturate.  Following  the 
suggestion  of  Hill  [1984],  Siscoe  et  ai .  [2002]  postulated 
that  a  saturation  potential  (FS  contributes  to  <FPC 


<f,pc  "  <C+X  “  ('  +§ij)  8($o  +  U)VswBTSin20 


=  $s)g 


(<Ko 


Q 

+  L<jVsw  BjStrr 


9 


(B3) 


Note  that  when  <J>E  C  <FS,  <t>PC  — *  4>E,  and  the  standard 
linear  relation  equation  (B2)  is  retrieved.  Sinee  the  ratio  4>E/ 
<E>s  >  0,  <f>s)  <  1 .  We  define  the  w  idth  of  the  storm  time 

geoeffeetive  gate  as  LGS  =  $s)  Sinee  eVs  =  $p<y 

2LYRE,  the  ratio  erfevs  %  2LY/LGS.  Beeause  LGS  is  a 
function  of  S\  the  ratio  of  interplanetary  to  magnetosphene 
eleetrie  fields  is  inherently  nonlinear.  No  single  ratio 
describes  all  situations. 

[55]  Siscoe  et  a l  [2002]  deriv  ed  an  expression  for  <J>S  in  k  V. 


<M  kV)  = 


1600  t/Psw(nPa) 
Y*p(mho) 


(B4) 


where  PSw  is  the  dynamte  pressure  of  the  solar  w  ind  and 
is  the  effective  Pedersen  conductance  of  the  polar  iono¬ 
sphere,  Oberet  ai.  [2003]  used  data  from  the  SWEPAM  and 
MFI  on  ACE  to  compare  <I>PC  predictions  of  the  S-H  and 
other  models  with  DMSP  measurements  during  the 
magnetic  storm  of  Mareh  2001.  Under  prestorm  conditions 
all  of  the  models  predicted  4>PC  values  measured  by  DMSP. 
During  the  main  phase  S-H  predictions  of  <I>PC  provided  an 
excellent,  upper  bound  envelope  for  DMSP  observations. 
Linear  models  predteted  mueh  larger  than  observed  <FPC. 
Burke  et  ai.  [2007b]  found  that  setting  £P  =  10  mho 
provided  excellent  agreement  between  calculated  and 
observed  <Fpc  during  a  large  number  of  magnetic  storms. 

[56]  In  computations  of  5Vs  we  used  hourly  averaged 
values  of  the  solar  wind  density  and  veloetty  from  SWE¬ 
PAM  and  IMF  components  in  solar-magnetospherie  coor¬ 
dinates  from  the  MFI  on  ACE.  In  applications  to  Siseoe-Hill 
we  set  4>o  =  25  kV,  LG  =  3.5  and  DP  =  10  mho. 

[57]  Acknowledgments.  Support  for  the  presented  work  was  provid¬ 
ed  by  the  Air  Force  Office  of  Scientific  Research,  Task  231 1SDA3,  and  AF 
contract  FA87I8-08-C-0012  with  Boston  College  The  authors  thank  Brian 
D.  Taplcy  of  the  University  of  Texas  for  use  of  the  reduced  GRACE 
measurements  that  form  the  empirical  basis  of  this  report. 

[58]  Amitava  Bhattachaijcc  thanks  the  reviewers  for  their  assistance  in 
evaluating  this  paper. 

References 

Bowman,  B  R.,  W.  K  Tobiska,  and  F.  A.  Marcos  (2006).  A  new  empirical 
thermospheric  density  model  JB2006  using  new  solar  indices,  AIAA 
Astrodvnamics  Conference .  Keystone,  Colo.,  AIAA  2006-6166. 

Bowman,  B  R.,  W.  K.  Tobiska,  F.  A.  Marcos,  C.  Y  Huang,  C.  S.  Lin,  and 
W.  J.  Burke  (2008),  A  new  empirical  thermospheric  density  model 


JB2008  using  new  solar  and  geomagnetic  indices,  Astrodynamic  Specia¬ 
list  Conference ,  Honolulu,  Hawaii,  AIAA  2008-6428. 

Bruinsma,  S.,  and  R.  Biancalc  (2003),  Total  densities  derived  from  accel¬ 
erometer  data,  J.  Spacecraft  Rockets ,  40,  230-236 

Bruinsma,  S.,  D.  Tamagnan,  and  R.  Biancalc  (2004),  Atmospheric  densities 
derived  from  Cl  IAMP/STAR  accelerometer  observations.  Planet.  Space 
Sci.,  52,  297-312. 

Burke,  W  J  (2007),  Penetration  electric  fields:  A  Volland-Stcm  approach, 
J.  Atmos  Sol  Terr.  Phvs.,  69,  1114-1126. 

Burke,  W.  J.  (2008),  Stormtime  energy  budgets  of  the  global  thermosphere, 
in  Mid-Latitude  Ionospheric  Dynamics  and  Disturbances,  Geophvs. 
Monogr.  Ser ,  vol.  181,  edited  by  P.  M.  Kintncr  et  al.,  pp.  235  246, 
AGU,  Washington,  D.  C. 

Burke,  W.  J.,  D.  R  Wcimcr,  and  N.  C.  Maynard  (1999),  Geoeffeetive 
interplanetary  scale  sizes  derived  from  regression  analysis  of  polar  cap 
potentials,  J.  Geophvs.  Res.,  104,  9989. 

Burke,  W.  J.,  C.  Y  Huang,  F.  A.  Marcos,  and  J  O.  Wise  (2007a),  Inter¬ 
planetary  control  of  thermospheric  densities  during  large  magnetic 
storms,  J.  Atmos.  Sol.  Terr.  Phvs.,  69,  279-  287. 

Burke,  W.  J.,  L.  C.  Gentile,  and  C.  Y  Huang  (2007b),  Penetration  electric 
fields  driving  main-phase  Dst,  J.  Geophvs.  Res.,  112,  A07208, 
doi:  1 0. 1 029/2006J  A0 12137. 

Burton,  R.  K  ,  R.  L.  MePhcrron,  and  C.  T  Russell  (1975),  An  Empirical 
relationship  between  interplanetary  conditions  and  Dst,  J.  Geophvs.  Res., 
80,  4204. 

Carovillano,  R.  L.,  and  G.  L.  Siscoe  (1973),  Energy  and  momentum  theo¬ 
rems  in  magnetosphene  processes.  Rev.  Geophvs.,  //,  289. 

Casali,  S.  J.,  and  W  N.  Barker  (2002),  Dynamic  calibration  atmosphere 
(DCA)  tool  for  the  high  accuracy  satellite  drag  model  (IIASDM),  AAS/ 
AIAA  Astrodvnomics  Specialist  Conference,  Monterey,  Calif.,  AIAA 
2002-4888. 

Ejiri,  M.  (1978),  Trajectory  traces  of  charged  particles  in  the  magneto¬ 
sphere,  J.  Geophvs.  Res.,  83,  4798. 

Hill.  T.  Wr.  (1984),  Magnetic  coupling  between  solar  wind  and  magneto¬ 
sphere:  Regulated  by  lonosphcne  conductance,  Eos  Trans.  AGU,  65, 
1047. 

Huang,  C.  Y,  and  W.  J.  Burke  (2004),  Transient  sheets  of  field-aligned 
current  observed  by  DMSP  during  the  main  phase  of  a  magnetic  storm, 
J  Geophvs.  Res.,  109,  A06303,  doi:  10. 1029 '2003J  A0 1 0067. 

Jacchia,  L.  G.  (1977),  Thermospheric  temperature,  density  and  composi¬ 
tion:  A  new  model,  SAO  Spec.  Rep.  375,  15  March. 

Kmpp,  D.  J ,  W.  K  Tobiska,  and  B  A.  F.mcry  (2005),  Direct  and  indirect 
thermospheric  heating  sources  for  solar  cycles  21  -23,  Sol.  Phvs.,  224, 
495  505. 

Marcos,  F.  A.,  M  J  Kendra,  J  M.  Griffin,  J.  N.  Bass,  D  R.  Larson,  and  J.  J 
Liu  (1998),  Precision  low  Earth  orbit  determination  using  atmospheric 
density  calibration,  J.  Astronaut.  Sci.,  6,  395  409. 

McComas,  D.  J.,  S.  J.  Bamc,  P  Barber.  W.  C.  Field  man,  J  L  Phillips,  and 
P.  Riley  (1998),  Solar  wind  electron,  proton,  and  alpha  monitor  (SWE¬ 
PAM)  on  the  Advanced  Composition  Explorer,  Space  Sci.  Rev.  86,  563 
612. 

Mlynezak,  M  G.,  ct  al.  (2005),  Energy  transport  in  the  thermosphere  during 
the  solar  storms  of  April  2002,  J.  Geophvs.  Res..  HO,  AI2S25, 
doi:  1 0. 1 029/2005JA0 1 1 1 4 1 . 

Obcr,  D.  M.,  N.  C.  Maynard,  and  W.  J.  Burke  (2003),  Testing  the  Hill 
model  of  transpolar  potential  saturation,  J.  Geoplns.  Res..  I08( A 12), 
1467,  dor  1 0. 1 029/2003 J  A0 10154. 

Roelof,  E.  C.,  and  D.  G.  Sibcck  ( 1993),  Magnetopause  shape  as  a  bivariate 
function  of  interplanetary  magnetic  field  Bz  and  solar  wind  dynamic 
pressure,  J.  Geophvs.  Res.,  98,  21,241 

Siscoe,  G.  L.,  G.  M.  Erickson,  B.  U.  6.  Sonncrup,  N.  C.  Maynard,  J.  A. 
Schocndorf,  K  D.  Sicbcrt,  D.  R.  Wcimcr,  Wr  W.  White,  and  G.  R.  Wilson 
(2002),  Hill  model  of  transpolar  potential  saturation  Comparison  with 
MUD  simulation,  J  Geophvs.  Res.,  I07( A6).  1075,  doi:  10. 1029/ 
2001 JA000I09. 

Smith,  C.  W  ,  M.  H.  Acuna,  L.  F.  Burlaga,  J  L.  LTIcurcux,  N.  F  Ness,  and 
J  Sehcifclc  (1998),  The  ACE  magnetic  field  experiment.  Space  Sci.  Rev., 
86,  613  630. 

Stem,  D  P  (2005),  A  historical  introduction  to  the  ring  current,  in  The 
Inner  Magnetosphere:  Physics  and  Modeling,  Geophvs.  Monogr.  Ser  . 
vol  155,  edited  by  T.  I  Pulkkincn,  N  A.  Tsygancnko.  and  R.  H.  Wr. 
Friedcl,  p.  I,  AGU.  Washington,  D.  C. 

Taplcy,  B  D.,  S.  Bcttadpur,  M.  Watkins,  and  C.  Rcigbcr  (2004),  The 
gravity  recovery  and  climate  experiment:  Mission  overview  and  early 
results,  Geophvs.  Res  Lett.,  31,  L09607,  doi:10.l029  '2004GL0 19920. 

Tcmcrin,  M.,  and  X.  Li  (2002),  A  new  model  for  the  prediction  of  Dst  on 
the  basis  of  the  solar  w  ind,  J.  Geophvs.  Res.,  I07( A 1 2),  1472, 
doi:  10.1029/2001  JA007530. 

Tcmcrin,  M.,  and  X.  Li  (2006),  Dst  model  for  1995  2002,  J  Geophvs. 
Res..  Ill,  A0422I,  doi:l0.l029'2005JA0l  1257. 


II  of  12 


A06306 


BURKE  ET  AL.:  STORM  TIME  THERMOSPHERE 


A06306 


Wcimcr,  D  R.  (2005),  Improved  ionospheric  clcctrodynamic  models  and 
application  to  calculating  Joule  heating  rates,  J.  Geophvs.  Res,,  110, 
A05306,  doi:10.!029/2005JA0l0884. 

Wcimcr,  D.  R.,  D.  M.  Ober,  N.  C  Maynard,  W  J  Burke,  M.  R.  Collier,  D.  J 
McComas,  T.  Nagai,  and  C.  W.  Smith  (2002),  Variable  time  delays  in  the 
propagation  of  the  interplanetary  magnetic  field,  J ,  Geophvs .  Res.. 
107( A8),  1211,  doi:  10.1 029/200 1 J A009 1 02. 

Wilson.  G.  R  ,  D  R,  Wcimcr,  J.  O.  Wise,  and  F.  A.  Marcos  (2006), 
Response  of  the  thermosphere  to  Joule  heating,  J.  Geophvs.  Res.,  Ill, 
A 1 03 1 4,  doi:  10. 1029/2 006J A0 1 1274 
Wright,  D  (2007),  Space  debris,  Phys.  Today ,  60(10),  35  40 
Wygant.  J.,  D.  Rowland,  H.  J.  Singer,  M.  Tcmcrin,  F.  Mozer,  and  M.  K. 
Hudson  (1998).  Experimental  evidence  on  the  role  of  the  large  spatial 
scale  electric  field  in  creating  the  ring  current,  J.  Geophvs.  Res ,  103, 
29,527. 


Yang,  Y.-H.,  J.  K.  Shoa.  A.  V.  Dmitriev,  C.-H.  Lin,  and  D.  M.  Ober  (2003), 
Saturation  of  IMF  Bz  influence  on  the  position  of  dayside  magnetopause, 
J  Geophys.  Res.,  108{ A3),  1104,  dor.10  1029  2002JA00962I. 

Young,  H.  D,  (1992),  Physics,  8th  cd.,  pp.  874  878,  Addison  Wesley, 
Reading,  Mass. 


W.  J.  Burke,  C  Y.  Huang,  C.  S.  Lin,  F.  A.  Marcos,  and  J.  O.  Wise,  Air 
Force  Research  Laboratory,  Space  Vehicles  Directorate.  29  Randolph  Road, 
Hanscom,  MA  01731-3010,  USA. 

L.  C.  Gentile  and  M  P,  Hagan,  Institute  for  Scientific  Research,  Boston 
College,  402  St.  Clement's  Hall,  140  Commonwealth  Avenue,  Chestnut 
Hill.  MA  02467-3862,  USA. 

D.  R  Wcimcr,  Virginia  Tech  National  Institute  of  Aerospace,  100 
Exploration  Way,  Hampton,  VA  23666.  USA. 


12  of  12 


